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Comparison  of  phenomenological  models  with  a microscopic  theory  for 
semiconductor  optical  nonlinearities 

I.  Rumyantsev t,  N.  H.  Kwongtt,  R.  Takayamaft  and  R.  Binderf 
t Optical  Sciences  Center,  University  of  Arizona,  Tucson,  AZ  85721  ,USA 
t Cooperative  Excitation  Project,  ERATO,  Japan  Science  and  Technology  Corporation 

Abstract.  A microscopic  theory  of  the  third  order  coherent  semiconductor  optical  response  is  used 
for  the  investigation  of  four-wave  mixing  signals.  The  theory  is  based  on  the  dynamics-controlled 
truncation  formalism  and  evaluated  for  the  case  of  resonant  excitation  of  heavy-hole  excitons.  This 
approach  allows  for  a direct  comparison  with  other  theories  (microscopic  or  phenomenological), 
and  yields  a detailed  understanding  of  important  many-body  effects  such  as  excitation-induced 
dephasing. 


Introduction 

The  nonlinear  optical  response  of  semiconductors  quantum  wells  has  long  been  an  impor- 
tant aspect  of  both  theoretical  and  experimental  semiconductor  physics  [1,  ].  The  un- 
derstanding of  many-body  effects  such  as  excitonic  effects  and  charge  carrier-correlations 
is  important  for  a correct  interpretation  of  nonlinear  optical  response  of  quantum  wells. 
Important  examples  of  recent  experiments  include  the  observation  of  biexcitonic  beatings 
in  four-wave  mixing  (FWM)  experiments  [ ] and  strong  polarization  dependence  of  the 
FWM  signal  in  semiconductor  microcavities  [ ]. 

The  theoretical  description  of  optically  excited  semiconductors  is  usually  based  on  var- 
ious modifications  and  extensions  of  the  so-called  semiconductor  Bloch  equations,  which 
are  equations  of  motion  for  the  wavevector-dependent  interband  polarization  function  and 
the  charge  carrier  (i.e.  electron  and  hole)  distribution  functions.  The  simplest  approx- 
imation is  the  Hartree-Fock  (HF)  approximation,  in  which  all  correlation  functions  are 
factorized  in  terms  of  one-particle  functions. 

However,  there  exist  several  theoretical  approaches  for  dealing  with  many-body  pro- 
cesses beyond  the  HF  approximation.  One  of  them  is  the  dynamics-controlled  truncation 
(DCT)  formalism,  which  is  a rigorous  approach  in,  for  example,  the  / (3)  nonlinear  regime 
if  one  restricts  the  analysis  to  the  coherent  limit  with  zero  electron  and  hole  density  initial 
conditions  | ].  These  equations  can  be  written  in  terms  of  the  interband  polarization  func- 
tion and  the  two-exciton  correlation  function.  If  consistent  with  the  excitation  conditions, 
the  equations  can  be  restricted  to  contain  only  optically  excited  ls-excitons.  In  that  case 
it  is  advantageous  to  expand  all  functions  in  terms  of  exciton  functions  and  keep  only  the 
smallest  set  of  excitons,  that  is  consistent  with  the  excitation  conditions  | ]. 

While  the  numerical  complications  only  allow  for  an  inclusion  of  a few  eigenstates  (Is 
in  our  case),  this  is  still  well  justified  approximation  for  a resonant  excitation,  especially 
in  semiconductor  with  large  exciton  binding  energies  are  used  or  spectrally  narrow  pulses 
(e.g.  picosecond  pulses).  From  this  approach  one  can  extract  and  calculate  the  nonlinear 
optical  response  function  and  extract  different  quantities  related  to  many-body  interaction 
effects,  e.g.  phase-space  filling,  Hartree-Fock  terms,  coherent  biexciton  contributions  and 
exciton  continuum  scattering  contributions.  These  quantities  can  be  directly  compared  with 
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the  phenomenological  parameters  entering  other  models  such  as  weakly  interacting  Boson 
(WIB)  model  (see  e.g.  [ ])  or  various  few-level  models.  The  comparison  with  the  W1B 
model  is  done  elsewhere  [ ],  while  the  comparison  with  few-level  models  is  discussed  in 
this  papers. 


1.  Theory  and  discussion 

In  this  section  we  discuss  both  the  microscopic  theory  based  on  the  DCT  formalism  as  well 
as  the  3rd-order  nonlinear  optical  response  based  on  a commonly  used  phenomenological 
five-level  model. 

The  five-level  system  corresponds  to  two  two-level  systems  representing  the  heavy  holes 
in  the  angular  momentum  states  [ ± 3/2  >,  the  electrons  in  the  spin  states  | ± 1/2  >,  and 
a phenomenologically  included  biexciton  state. 

Using  an  atomic  density  matrix  approach  one  can  write  the  general  equation  of  motion 
for  the  density  operator  in  the  form 

ifi  — Pmn  = ( hO->mn  ( Y ^ Pi i ) Pmn  T ^ ' ( E ju  Pirn  PniHjm),  (1) 

i>g  i 


where  the  interaction  hamiltonian  H'mn  = — • (E  + LP),  fimn  is  the  dipole  matrix 
element,  E is  the  optical  field,  and  P is  the  spin-resolved  polarization  function.  The  many- 
body  effects  are  accounted  for  by  the  inclusion  of  the  biexciton  level,  the  excitation-induced 
dephasing  (EID)  associated  with  the  parameter  y',  and  the  local  field  factor  L.  Iterating 
this  equation  to  the  third  order  one  can  derive  the  expression  for  the  total  polarization  using 
p = —eT r(rp)  = J2nm  P-nmPmn-  For  the  spin  +1  polarization  component  one  obtains 
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A similar  expression  can  be  derived  for  the  spin  —1  polarization  component.  We  use  the 
notation  p±  to  denote  the  first  order  density  matrix  (linear  optical  response)  E±  = jj,±E, 
L = L\h\2/V,  where  V is  the  volume  or,  as  in  our  case  of  a quasi-two  dimensional  system, 
the  total  area  of  the  system. 
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Fig.  1.  Four-wave  mixing  signal  vs  pump-probe  delay  time  in  co-polarized  (dashed  lines)  and 
cross-polarized  (solid  line)  configuration  computed  with  the  microscopic  (a)  and  phenomenoligical 
(b)  theory. 


We  now  turn  to  the  microscopic  formulation.  The  equation  for  the  third-order  polariza- 
tion that  one  can  be  obtain  from  the  DCT  equations  in  terms  of  exciton  wave  functions  can 
be  written  in  the  following  form: 


pfit) 
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Here,  p is  the  first  order  coherence  function.  The  two  terms  on  the  first  line  of  this  equation 
represent  the  phase  space  filling  (PSF)  and  Hartree-Fock  (HF)  contributions.  This  is 
followed  by  the  bound  biexciton  term  on  the  2nd  line,  (restricted  to  the  singlet  channel) 
and  the  rest  of  the  equation  contains  the  continuum  scattering  terms,  with  contributions  in 
both  triplet  and  singlet  channels.  Here  “singlet”  (/o  and/-)  and  “triplet”  (/+)  refer  to  the 
electron  spin  states.  The  definitions  of  the  parameters  C]  v and  Vst  (see,  e.g.,  | ]),  as  well 

as  /o,  and  f~(t  — t")  are  lengthy  and  not  given  here,  they  include  multiple  integrals  with 
exciton  eigenfunctions,  the  Coulomb  potential.  The  biexciton  contribution  and  continuum 
contribution  in  the  singlet  channel  contain  also  the  biexciton  wave  function  (which  can  be 
found  in  Ref.  [ ]). 

Equations  (3)  and  (2)  are  a suitable  basis  for  a detailed  comparison  of  parameter  enter- 
ing the  phenomenological  model  with  the  interaction  terms  derived  from  the  microscopic 
theory.  Whereas  a strict  one-to-one  correspondence  does  not  exist,  one  can  indeed  identify 
approximate  correspondences  and,  thus,  validate  to  a certain  extent  the  phenomenolog- 
ical model.  As  far  as  numerical  agreement  is  concerned,  we  show  in  Fig.  1 results  for 
time-integrated  four-wave  mixing  signal  in  a ZnSe  semiconductor  quantum  well. 

In  conclusion,  we  present  a theoretical  analysis  of  ultrafast  optical  nonlinearties  in 
semiconductors.  We  believe  that  our  detailed  comparison  of  phenomenological  model 
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with  a microscopic  theory  is  important  in  order  to  validate  conclusions  about  physical 
processes  obtained  from  the  widely  used  phenomenological  models. 
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